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Abstract 

The three- spin chain with Heisenberg XY- interaction is simulated in a three- qubit nuclear 
magnètic resonance (NMR) quantum computer. The evolution caused by the XY- interaction 
is decomposed into a series of single- spin rotations and the J- coupling evolutions between the 
neighboring spins. The perfect state transfer (PST) algorithm proposed by M. Christandl et al 
[Phys. Rev. Lett, 92, 187902(2004)] is realized in the XY- chain. 
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I. INTRODUCTION 



In 1982, R. P. Feynman proposed the idea of quantum computer, and pointed out that 
a quantum computer can simulate physics more efficiently than its classical counterpart Q . 
The information carriers for quantum computation are qubits. Unlike a classical bit, a qubit 
can lie in a superposition of two states, according to the quantum mechanical principle of 
superposition. Essentially, the superpositions of quantum states lead to the advantages of 
the quantum computers over the classical computers. The quantum gates and quantum 
networks proposed by D. Deutsch provi de a convenient method for people to think about 
how to build a quantum computer in a similar way to build a classical computer 

as 

. His 



work can be thought as a milestone in the history of quantum computer 



The current 



quantum network theory has shown that it is possible to construct an arbitrary n-qubit 



quantum gate by using only a finite set of single- qubit gates and two- qubit gates 



that these bàsic quantum gates are universal for quantum computation 



qubit ; 

IBI 



fll 



so 



The prime 



factorization algorithm proposed by P. W. Shor and the quantum search algorithm 
proposed by L. K. Grover show the potential advantages of the quantum computers and 
accelerate the development of quantum computation. 

There are several physical systems that can implement quantum computation. They 



are liquid nuclear magnètic resonance (NMR) 11211. quantum dots 



nqi 
[i 



16l . Il7j . electron spins [TJ 



& 



13, 



Q, 



so 



l20j . superconduction qubits 



id NMR 
21, U, 



and 



trapped ions 

cavity QED systems [22|, |24j, |25j, [2a| . Because of its technologic maturation and convenience 
in manipulation, liquid NMR has been an important experimental method to implement 



quantum algorithm, error-correcting code, and simulate quantum systems 
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In the above systems, the interactions between qubits, at least between the neighboring 
qubits, are necessary for quantum computation. The Heisenberg interaction naturally exists 
in the various spin systems. In the liquid NMR system, the Heisenberg interaction exists 
in form of Ising interaction j^J. In the other systems, the Heisenberg interaction takes 
more various forms. D. P. DiVincenzo et al pointed out that the Heisenberg interaction 
alone can be universal for quantum computation, if the coded qubit states are introduce 
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This result is exciting, because the single- spin operations, which usually cause 



additional difficulties in manipulations in some systems, can be avoided. The perfect state 
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transfer (PST) algorithm proposed by M. Christandl et al satisfies such a condition that 



no single-spin operations are needed 



44 1 . The algorithm can transfer an arbitrary quantum 



state between the two ends of a spin chain or a more complex spin network in a fixed period 
time only using the XY- interaction. If the state is transferred in a more than three spin 
chain, the coded qubits are needed, so that the chain is extended to a network. Compared 
with the state transfer based on SWAP operations, where single-spin operations are used to 
switch on or off the couplings between spins [45], the PST algorithm is easy to implement 
in some sòlid systems. 

The Heisenberg interaction is expected to play an important role in building large- scale 
quantum computers, and it has become an interesting tòpic in the field of quantum informa- 
tion. M. C. Arnesen et al' work indicated that t he q uantum entanglement phase transition 
occurs in the one dimension Heisenberg model [46]. L. Zhou et al pointed out that the 
thermal entanglement can be enhanced in an anisotropic Heisenberg XYZ chain 47|- J. P. 
Keating et al separated a quantum spin chain into two parts, and computed the entropy 
of entanglement between them j3]. M. Mohseni et al proposed a fault-tolerant quantum 
computation using Heisenberg interactions J^. The other issues related to the Heisenberg 
model, such as the Heisenberg chain with the next-nearest-neighbor interaction, are also 



discussed 
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531 ] . In experiment, the quantum entanglement phase transition in a 
two- spin Ising- chain was demonstrated in an NMR quantum computer j^. 

)uantum simulation has been an interesting tòpic since the quantum computer is bornQ, 



Quantur 



j. Liquid NMR has displayed its powerful ability to simulate quantum systems, and 
various of quantum systems have been successfully simulated in NMR quantum computers 
In this paper, we simulate a three- spin XY- chain with the Heisenberg 
interaction and realize the perfect state transfer(PST) algorithm 44j using a liquid NMR 
quantum computer. 



II. SIMULATING THE THREE- SPIN XY- CHAIN USING LIQUID NMR 

The Hamiltonian for a three spin XY- chain with the neighboring Heisenberg interaction 

is 



hxy = ^j(°y x + °y y + k4 + °y v 



(i) 
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where & x / y {j = 1,2,3) are the Pauli matrices for the angular momentum of the spins, and 
J is the coupling constant between two spins. For convenience in expression, h has been set 
to 1. The evolution caused by Hxy can be expressed as 

U(t) = e~ iHxY \ (2) 

where t is the evolution time. In order to represent U(t) as a liquid NMR version, we 
introduce two commutable operators A = (cr^cr| + <7^<7^)/2, and B = (cTyCTy + a x a x )/2. U(t) 
can be rewrite as U(t) = U A (t)UB(t), where 

U A (t) = e~ iJtA = e -^l°l+°>Dl\ (3) 

U B (t) = e ~ iJtB = e - ijt (<rl°i+°i°i)/2 m (4) 

We define three operators L x = a\a 2 x /2, Ly = a^a^/2, and L A = o 1 x o\o y \/2. These three 
operators can be viewed as the three components of the angular momentum vector denoted 
by L A , because they satisfy the commuting conditions [L x , Ly} = iL A , [Ly, Lf] = iL x , and 
[Lf, L£] = iLy. Eq. Q can be rewritten as 

U A (t) = e - utiL £ +L ^ = e -*(^)L A - (5) 

where the vector n = (l/y/2, \/\[2, 0), and it denotes the direction of the rotation axis for 
{7,4 (i). The separate angles between n and x, y, z axes are n/4, n/4, and n/2, respectively. 
Using the theories of angular momentum, we obtain 



U A (t) = e -iïLÍ e -iV2JtL£ e i*LÍ 
Jt s 1 , Jt 



In a similar way, through defining L B = a 2 a x /2, Ly = <j^<j 2 /2, and Lf = a^alal/2 as 
the three components of the angular momentum vector denoted as L B , we obtain 



U B (t) = e -^ e -*V2JtLÏ^L? 
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Jt , i . / Jt \ s x 2 i 2 :■! 



COS 



72 )_ 72 sin( 72 )((Tï/ ^ + ÍT:E(Ta;) ' (7) 



One can prové the last equations in Eqs.flHJ) and Q directly through Eqs. (JHJ) and (JIJ) using 

(((7^2 + ^3)/^!! = 1; and ((^2 + ^3)^)2 = L Combining Eqg . Q and W e 

obtain 



U(t) = e^^^e^^'^e^^e^^^e^^^e^^. (8) 

Each of the six factors in Eq. (jHJ) can be realized using liquid NMR. Consequently the three- 
spin XY- chain can be simulated in a three- spin liquid NMR system. 



III. IMPLEMENTING THE PERFECT STATE TRANSFER ALGORITHM IN 
THE XY- CHAIN 

The PST algorithm was proposed by M. Christandl et al 44], and it can be implemented 
in the XY chain. The algorithm can transfer an arbitrary quantum state between the two 
ends of the chain in a fixed period time, only using the XY- interaction. Unlike the state 
transfer based on SWAP operations [45], the PST algorithm do not require single- spin 
operations. Hence the algorithm is more feasible to realize in some systems, such as the 
electron-spin-resonance system, where the single- spin operations cause many experimental 
difficulties Q 

Letting (p = Jí/\/2, Eq. (JHJ) is represented as the matrix 



U(t) = 



(l 

















o o N 





COS 2 if 


-^sin(2^) 





— sin 2 <p 











-^sin(2y?) 
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— sin 2 <p 


-^sin(2y?) 





cos 2 ip 

















-^sin(2^) 





cos(2<£>) 


-^sin(2^) 











— sin 2 tp 





-^sin(2^) 


cos 2 ip> 


1° 

















1 



(9) 
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The order ofthe basis states is |000), |001), |010), |011), |100), |101), |110), |111), where |0) 
and |1) denote the spin up and down states, respectively. When t = -^j, one obtains 



7Ï 



V2J J 









0^ 




0-10 










0-10 













-1 







-10 










0-1 










0-10 




















|100), Z7|010) = -|010), 


Ï7|011) = -|110) 


|011), and U\lll) = 


1111). 


We use |V>)in 



(10) 



(a|0) + 



V 

Obviously, C/|000) = |000), Z7|001) = 
-|001), C/|101) = — 1101>, C7|110> = 
/3|1))|00) as the input state by setting spin 1 into state (a|0) + /3|1)), where a, (3 are two 
arbitrary complex numbers. U(-j=j) transforms \ip)i n to |00)(a|0) — /3|1)), where spin 3 
lies in state (a|0) — and the perfect state transfer is completed. Obviously, a simple 

operation a z can transform (a|0) — (3\1)) to («|0) + /3|1)), and therefore \ip)i n is transferred 
from spin 1 to spin 3. 

The implementat ion of PST algorithm in two- or three- spin chain does not require the 
coded qubits. However in more than three-spin chain, the coded qubits are needed to design 
so as to extend the chain to a more complex network. The details can been found in 



IV. REALIZATION IN A THREE- NMR QUANTUM COMPUTER 



The experiment uses a sample of Carbon-13 labelled trichloroethylene (TCE) dissolved in 
d-chloroform. Data are taken with a Bruker DRX 500 MHz spectrometer. The temperature 
is controlled at 22°C. 1 H is denoted as qubit 2, the 13 C directly connecting to 1 H is denoted 
as qubit 1, and the other 13 C is denoted as qubit 3. The three qubits are denoted as Cl, H2 
and C3. The Hamiltonian of the three-qubit system is 



H 



xmr = --kv x o\ - iïv 2 o 2 z - nu 3 a 3 z + -ix J 12 (J l z (J 2 z + -^^23^^ + -7rJ 13 <7^af, 



11^ 
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where ui, z/2, z/3 are the resonance freqüències of Cl, H2 and C3, and z/3 = v\ + 904.4Hz. The 
coupling constants are measured to be J12 = 200.9 Hz, J 2 3 = 9.16 Hz, and J13 = 103.1Hz. 
The coupled-spin evolution between two spins is denoted as 

[ Tjl ] = e- i ** J i tTai ° X ; (12) 

where l = 1,2, 3, and j ^ l. [tji] can be realized by averaging the coupling constants other 
than Jji to zero p 71] . 

The three- body and two- body interactions in Eq. © can be expressed as 

í-íxW 3 -i-a 1 i-a 3 i-a 1 a 2 a 3 i-a 1 -i-a 3 1 1 a\ 

^ L s x 2 y = g 4 v g 4 ^ g 8 2 2 U 2 ç 1 4 yg (14) 

g l 8 a y' T z cr x — e é xe * »e % % a * tr * a *e 7 x e * y , (15) 



e *f^J^o2 — e *i^ e -íi ^gíf^^^g-if^gíf^^ (is) 



g-^Sg-^g-^í^V^e 1 * ». (18) 
Through substituting Eqs. (jl3H18|) into Eq. (jSJ), and after simplification, one obtains 



xe *f (^+<r|) e -if CTjg-if <ï\a\al e -i\ a% e ~i<pa 2 z a 3 z ^ ^a\al<j% e -i\a\ e %\a% ^ 

In Eq. (fTÜ|) . e l ï^ is realized by a 7r/2 radio frequency (rf) pulse exciting H2 along y-axis. 
Such a pulse is denoted by The operation g*f ( <T ^+ f7 S) j s realized by a nonselective 

pulse [7r/2]^' 3 , exciting Cl and C3 simultaneously The widths of [7i"/2]^ and [7r/2]^' 3 are so 
short that they can be ignored. e 1 ^ ( <T z +<7 *) is realized by a pulse sequence 
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K' 3 - H' 3 , (20) 

where the time order is from left to right. The operation selective for Cl or C3 can be 
realized by the established pulse sequence For example, e 1 * ^ is realized by 



[\\l = [\\Y-^° l -[-\\Y- (21) 



According to C. H. Tseng et al's work |60|, e~n< a * a * is realized by 



7T,n , ,o r 9 , , IX , r, r 1 , r 7T . r 9 , ,71 



\ — — \ Z - -7T 



] " M - [tt] - ío) - - (22) 



and e^ * * * is realized by 

hp - m; - [^i - ifi; - [^i - if is - i^i - [fiï. (23) 

The case of only one proton in the sample makes the three- body interactions realize much 
easy One should note that the direct coupling between Cl and C3 is not used. 
We choose the state 

PiniA = o\ (24) 
as the initial state to simulate the XY- chain in the three-spin system. The pulse sequence 



- [fií - \9rad\ M - (25) 
transforms the system from the equilibrium 

Pe Q = lc(l] + II) + 1hI 2 z , (26) 

to piniA jüï|' where jc an d Jh denote the gyromagnetic ratios of 13 C and 1 H, and [grad] z 
denotes a gradient pulse along z- axis. The irrelative overall factors have been ignored. 
Using [U B , PiniA] = 0, we obtain p A (t) = U (t) p iniA U ] (t) = U A (t)p iniA U A (t). In experiments, 
we replace U(t) by U A (t), in order to simplify experimental procedure, and obtain 



S 



p A {t) = a\ cos 2 (p + °Wx^ sin(2^) - o^oj sin 2 (27) 

When í = one obtains P^(-^j) = — <T l cr z <T yj which means that the state a y has been 
transferred from Cl to C3. Similarly, if the initial state is chosen as 

PiniB — V x i (28) 

we obtain 



p B (t) = U{t) PimB U\t) = U B {t)p miB U B {t) 

= o\ cos 2 <p + a l a l-^ sm(2<f) - o\o\ol sin 2 <p. (29) 

Obviously, Pb(^j) = —°'l 'ï 'l, which means that a x has been transferred from Cl to C3. 

We represent the results of the implementation by NMR spectra. When (p changes, the 
amplitudes of Cl and C3 change as cos 2 y9 and sin 2 yj, respectively. When the initial state 
is chosen as PmiA, the experimental results are shown as Fig. Q The data for Cl are 
marked by "+", and are fitted as A\ cos 2 tp; the data for C3 are marked by " * ", and are 
fitted as A3sm 2 (p. The two constants A\ = 6.20 and A% = 5.65, with arbitrary units. The 
experimental results, barring two data for Cl, show a good agreement with the theoretical 
expectations. Figs. El show the spectra when the state transfers occur. When ip = 0, 
(p = 7r/2, (p — 7r, <p — 37r/2, and p = 2n, the system lies in cr* (the initial state), — cr^cr 2 ^, 
<jy, —<Tl<T^(jy, and <Ty, shown as Figs. El^a-e), respectively. The experimental results, barring 
the signals of Cl in Figs. I2fb) and (d) of which amplitudes are shown in Fig^ agree with 
the theoretical expectation quite well. Theoretically, the signals of Cl in Figs. Efb) and (d) 
should not appear. The time duration for implementing Ua is about 200ms, which is in the 
same order with the decoherence time. Hence the decoherence time limit results in main 
errors. Moreover, the imperfection of the pulses and the inhomogeneity in the magnètic field 
also cause errors. The similar results can be obtained when the initial state is chosen as 
PiniB- Figs. El show the implementation of the perfect state transfer when the initial state 
is PiniB- When <p = and <p = n/2, the system lies in al (the initial state) and — <T*cr 2 <7^, 
shown as Figs. Efa-b), respectively. 
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V. CONCLUSION 



We have simulated the three- spin XY chain using liquid NMR. Through defining proper 
operators, we use the theories of angular momentum to decompose the evolution caused by 
XY- coupling into a series of factors that can be realized by rf pulses and J- couplings. Such 
an analogue can be helpful for solving the general problems on the Heisenberg chain. As 
an example for the application of the XY- chain in quantum computation, the perfect state 
transfer algorithm is realized in the chain. 

The evolution caused by XY- couplings can be represented by single- spin operations and 
the J- couplings, although there are no real XY- couplings in liquid NMR. In the sample used 
in our experiments, the coupling constants are not equal to each other. However we simulate 
the equal couplings in the XY- chain through choosing the proper evolution time. For the 
PST in more than three spin networks, the coupling strengths are needed to be designed 



m a proper manner 



44j . Our work has shown that such couplings are easy to simulate in 



NMR. All these facts represent the powerful function of the liquid NMR in implementing 
quantum computation. 
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FIG. 1: The graph of the amplitudes of Cl and C3 vs. ep = Jt/y/2. The amplitudes have arbitrary 
units. The data for Cl are marked by "+", and are fitted as A\ cos 2 ip; the data for C3 are marked 
by " * ", and are fitted as ^3 sin 2 ip, where A\ = 6.20 and A% = 5.65. 
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FIG. 2: The NMR spectra for the implementation of the perfect state transfer when the initial 
state is &y. When tp = 0, tp = vr/2, p = n, p = 3tt/2, and p = 2ir, the system lies in a y (the initial 
state), — a\a 2 z ay, a y , —aja^ay, and a y , shown as Figs. (a-e), respectively. Fig. (a) is the reference 
spectrum used to calibrate the phases of the signals in Figs. (b-e). 



15 




FIG. 3: The NMR spectra for the implementation of the perfect state transfer when the initial 
state is a\. When ip = and 92 = 7r/2, the system lies in o\ (the initial state) and —G\a 2 z a\ shown 
as Figs. (a-b), respectively. Fig. (a) is the reference spectrum. There is a n/2 phase difference 
between the signals in Fig. (a) and Fig. Efa). 
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